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_ We calculate the nuclear cross section for coherent and incoherent vector meson production within 

the QCD color dipole picture, including saturation effects. Theoretical estimates for scattering on 
both light and heavy nuclei are given over a wide range of energy. 
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l ■ I. INTRODUCTION 

One of the main theoretical expectations for the high energy regime of Quantum Chromodynamics (QCD) is the 
saturation of the parton densities in hadrons and nuclei at small values of the Bjorken-x variable and the formation 
of a Color Glass Condensate (CGC). This is one of the main topics of hadron physics to be explored in the new 
accelerators, such as the LHC and possibly the future electron-ion collider. The search for signatures of the CGC 
has been subject of an active research (for recent reviews see, e.g. \m\ with the experimental data being successfully 
described by phenomenological models based on saturation physics 0, y, 0, B S 0| • However, more definite conclusions 
are not yet possible. In order to discriminate between these different approaches and test the CGC physics, it would 
be very important to consider an alternative search. To this purpose, the future electron-ion colliders offer a promising 
opportunity [!,[§]. 

In a series of papers [10, LkH Qjl, [l3| w e investigated the prospects of observing the CGC in a future electron-ion 
Q\ [ collider (For related studies see Refs. [l4 ITU. [l6jV As it has been already emphasized in these papers, the advantage 
■ of using nuclear targets is that the saturation scale Q s is much larger and this is crucial for the observation of most of 
the CGC effects. After these studies our conclusion was that it is very difficult to disentangle CGC effects from the 
standard linear QCD looking only at the nuclear inclusive observables, such as the nuclear structure functions F^, 
F% ' A and Ffi. On the other hand, the study of diffractive observables was shown to be promising, as demonstrated 
in Refs. [Tll.Tl3j. In particular, in [l3| we found out that the saturation models are able to describe the current data 
on the nuclear structure function FA and predict that the contribution of diffractive events to the total cross section 
should be of ks 20% at large A and small Q 2 . In the asymptotic limit of very high energies diffractive events are 
expected to form half of the total cross section, with the other half being formed by all inelastic processes [13] • This 
observation motivates a more detailed study of diffractive processes. In particular, we can expect that the study of 
exclusive processes, such as diffractive vector meson production, can be useful to determine the QCD dynamics at high 
energies. It is important to emphasize that the experimental HERA data on vector mes on p roduction in ep processes 
are successfully described by the phenomenological saturation models (See, e.g., 0, [H, 13 US HH ) • ^ n contrast, 
very little is known about vector meson production off nuclei at high energies. Exclusive vector meson production 
in eA interactions can be classified as coherent or incoherent. If the reaction leaves the target intact, the process 
is usually called coherent. Otherwise it is called incoherent. These two types of diffractive hadron production were 
studied recently in Refs. [23|, [24|, HH, [27| using the CGC formalism. The main conclusion of these works is that 
coherent production can serve as a sensitive probe of the high energy dynamics of nuclear matter, while incoherent 
production measures the fluctuations of the nuclear color fields. These findings are an additional motivation for our 
study of nuclear vector meson production. We focus our analysis on exclusive vector meson production by real and 
virtual photons, which may help us to understand several physical issues, besides saturation effects. Among these we 
find, for example, the transition from the soft dynamics (at low virtualities of the photon Q 2 ) to the hard perturbative 
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regime at high Q 2 and the relative weights of coherent and incoherent interactions between the projectile and the 
target nucleus. 

In this paper we calculate the coherent and incoherent production cross sections of vector meson production using 
the dipole approach and a nuclear saturation model based on CGC physics. The main input for our calculation is the 
dipole-target cross section, cr^ gct = <7di P (x,r), which is determined by the QCD dynamics at small x. In the eikonal 
approximation it is given by: 

<Tdip{x, r) = 2 J d 2 bM(x, r, b) (1) 

where Af(x,r,b) is the forward dipole-target scattering amplitude for a dipole with size r and impact parameter b 
which encodes all the information about the hadronic scattering, and thus about the non-linear and quantum effects 
in the hadron wave function (see e.g. [l|). It can be obtained by solving the BK (JIMWLK) evolution equation in the 
rapidity Y = ln(l/a;). Many groups have studied the numerical solution of the BK equation, but several improvements 
are still necessary before using the solution in the calculation of observables. In particular, one needs to include the 
next-to-leading order corrections into the evolution equation and perform a global analysis of all small x data. It is a 
program in progress (for recent results see [28l. [29j). In the meantime it is necessary to use phenomenological models 
for Af which capture the most essential properties of the solution. Following [l3| we will use in our calculations the 
model proposed in Ref. [3(3], which describes the current scarce experimental data on the nuclear structure function 
as well as includes the impact parameter dependence in the dipole nucleus cross section. In this model the forward 
dipole-nucleus amplitude is given by 



Af A (x,r,b) = 1 - exp 



-^a dp (x,r 2 )T A (b) 



(2) 



where Od v is the dipole-proton cross section and T A (b) is the nuclear profile function, which is obtained from a 3- 
parameter Fermi distribution for the nuclear density normalized to A (for details see, e.g., Ref. [Uj]). The above 
equation, based on the Glauber-Gribov formalism [321 ] . sums up all the multiple elastic rescattering diagrams of the 
qq pair and is justified for large coherence length, where the transverse separation r of partons in the multiparton 
Fock state of the photon becomes a conserved quantity, i. e. the size of the pair r becomes eigenvalue of the scattering 
matrix. It is important to emphasize that for very small values of x, other diagrams beyond the multiple Pomeron 
exchange considered here should contribute (e.g. Pomeron loops) and a more general approach for the high density 
(saturation) regime must be considered. However, we believe that this approach allows us to estimate the magnitude 
of the high density effects in the kinematical range of the future eA colliders. 

In the present work, we shall use two models for the dipole-proton cross section a dp- One is the very popular 
GBW model Q, which interpolates between the small and large dipole configurations, providing color transparency 
behavior, a dp ~ r 2 , as r — > and constant behavior, a dp ~ o"0j a t large dipole separation. This model is no longer 
able to describe the most recent HERA data and it has been replaced by other parameterizations. However we shall 
keep using it as a baseline, which will be compared with other dipole-proton cross sections, giving us an estimate of 
the sensitivity of the observable to changes in the dipole cross sections. The parameterization of the dipole-proton 
cross section is given by the eikonal-like form [3j, 



Q s 2 (i)=(|) GeV 2 , (3) 

where the saturation scale Q 2 defines the onset of the saturation phenomenon, which depends on energy, and x = (Q 2 + 
4m 2 )/(Q 2 + W 2 N ) (See [3j for details). One of the drawbacks of the GBW model is that it has no impact parameter 
dependence. It is assumed that the impact parameter dependence of TV can be factorized as J\f(x, r, b) = Af(x, r)S(b) 
and this last function is integrated over 6, giving rise to the parameter <j . 

During the last years an intense activity in the area resulted in more sophisticated dipole proton cross sections, 
which had more theoretical constraints and which were able to give a better description of the more recent HERA 
and/or RHIC data [H HH, H, E| H H [p . In what follows we will use the b-CGC model proposed in Ref. [H, 
which improves the IIM model pi] with the inclusion of the impact parameter dependence in the dipole proton cross 
sections. The parameters of this model were recently fitted to describe the current HERA data [22j ■ Following [2(| 
we have that the dipole-proton cross section is given by: 

< GC (*,r 2 )^ / d 2 b^ (4) 



<4 BW (£, r 2 ) 



ctq 



1 — exp 



(x) r z 
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where 



da, 



dp 



d 2 b 



2N p (x,r,b) = 2 x 



2(7,+ 



ln(2/rQ s p 



1 _ exp _a ( br Qs) 



rQ s < 2 
rg s > 2 



(5) 



with y = ln(l/a;) and n = x"{ls)/x'(ls), where \ is the LO BFKL characteristic function. The coefficients a and 
b are determined uniquely from the condition that M p (x,r) and its derivative with respect to rQ s are continuous at 
rQ s = 2. In this model, the proton saturation scale Q s now depends on the impact parameter: 



exp 



2£?CGC / 



(6) 



The parameter -Bcgc was adjusted to give a good description of the i-dependence of exclusive J/ip photoproduction. 
Moreover the factors A/"o and 7 S were taken to be free. In this way a very good description of F 2 data was obtained. 
The parameter set which is going to be used here is the one presented in the second line of Table II of (22J: 7^ = 0.46, 
Bcgc = 7.5 GeV~ 2 , j\f = 0.558, x = 1.84 x 10~ 6 and A = 0.119. 

The paper is organized as follows. In next section (Section [TTJ) we discuss the coherent and incoherent vector 
meson production and present our main formulas. In Section IlIII we present our predictions for the energy and Q 2 
dependences of the p and J/^f total cross sections considering the coherent and incoherent cases for two different 
nuclei. Moreover, we estimate the ratio between the coherent and incoherent cross sections. Finally, in Section HVl we 
summarize our main conclusions. 



II. EXCLUSIVE VECTOR MESON PRODUCTION IN THE COLOR DIPOLE APPROACH 



In the color dipole approach the exclusive vector meson production (7*^ — > VY) in electron-nucleus interactions 
at high energies can be factorized in terms of the fluctuation of the virtual photon into a qq color dipole, the dipole- 
nucleus scattering by a color singlet exchange and the recombination into the final state vector meson. This process 
is characterized by a rapidity gap in the final state. If the nucleus scatters elastically, Y = A, the process is called 
coherent production. On the other hand, if the nucleus scatters inelastically, i. e. breaks up (Y = X), the process 
is denoted incoherent production. In this paper we will consider the color dipole description of the 7* A — > VY 
(V = p, J/ty) process, which is quite successful for the proton case [1, [H, [ill I20I I21I HI], [33] and can be extended 
to nuclei targets with the Glauber-Gribov formalism (For recent reviews see |38l . [39|). To a large extent we follow 
the pioneering t 
developments [4.' 
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where these issues were first addressed and also their further 
To be more precise, in this work we will update the calculations presented in [52] |. 
extending them also to electroproduction and treating more carefully the two production mechanisms, coherent and 
incoherent. For this we shall use the formalism described in detail in Ref. (49[. For the dipole cross section we shall 
use the parameterization introduced in [22| , which was able to give a very good description of data on exclusive vector 
meson production on proton targets. In this sense our work is an extension of [22j to nuclear targets. 

As discussed above, the exclusive vector meson production in electron-nucleus interactions can be classified as 
coherent or incoherent. If the reaction leaves the target intact, the process is usually called coherent, and the mesons 
produced at different longitudinal coordinates and impact parameters add up coherently. The corresponding integrated 
cross section is given in the high energy regime (large coherence length: l c ^> Ra) by [49( 



oh ( 7 *A 



VA) = J d 2 b|| j d 2 r J dz^>* v (r,z)Af A (x,r,b) * n 



(r,z,Q 2 )\ : 



(7) 



On the other hand, the diffractive incoherent production of vector mesons off nuclei, j*A — * VX, is associated with 
the breakup of the nucleus. In this case one sums over all final states of the target nucleus, except those that contain 
particle production. The t slope is the same as in the case of a nucleon target. Therefore we have: 



( 7 *A -» VX) 



\XmA(s, t = 0)\ 2 



I671-.B 



v 



where at high energies (l c Ra) [491 ]: 

\lmA(s, t = 0)| 2 = J d 2 bT A (b) I J d 2 r J dz ^* v {r, z) a dp exp[-i a dp T A (b)\ (r, z, Q 2 )' 2 



(8) 



(9) 
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V(m v ) 


t-i9 

m f R A 


N L N T if 


(MeV) 


(GeV) (GeV -2 ) 




p(776) 


0.14 12.9 


0.853 0.911 l/y/2 


J/* (3097) 


1.4 2.3 


0.575 0.578 2/3 



TABLE I: Parameters and normalization of the boosted-Gaussian overlap function. 



The qq pair attenuates with a constant absorption cross section, as in the Glauber model, except that the whole 
exponential is averaged rather than just the cross section in the exponent. The coherent and incoherent cross sections 
depend differently on t. As discussed in j43[, at small-i (— tR\/3 <C 1) coherent production dominates the leptopro- 
duction of vector mesons, with the signature being a sharp forward diffraction peak. On the other hand, incoherent 
production will dominate the vector meson production at large-i (—tR\/3 3> 1), with the ^-dependence being to a 
good accuracy the same as in the production off free nucleons. Therefore, it is expected that this signature will allow 
us to separate the two contributions even at a limited resolution in t. However, at an electron-ion collider, if the 
nucleus remains intact in a small-i interaction, it will escape too close to the beam to be detectable. This implies that 
the experimental separation between coherent and incoherent production will be a challenging task. In this paper 
we focus our study on the energy and Q 2 dependence of the total coherent and incoherent cross sections and their 
relative weights. We postpone the study of the i-dependence of these cross sections for a future publication. 

In the Eqs. ([7]) and (J9j) the functions Az, r) and 4>^(z, r) are the light-cone wavefunctions of the photon and 

vector meson, respectively. The quark and antiquark helicities are labeled by h and h and reference to the meson 
and photon helicities are implicitly understood. The variable r defines the relative transverse separation of the pair 
(dipole) and z (1 — z) is the longitudinal momentum fraction of the quark (antiquark). In the dipole formalism, the 
light-cone wavefunctions 4^ Az, r) in the mixed representation (r, z) are obtained through two dimensional Fourier 
transform of the momentum space light-cone wavefunctions 4^ A\z, k) (see more details, e.g. in Ref. [Hj]). The 
normalized light-cone wavefunctions of longitudinally (L) and transversely (T) polarized photons are given by [20I ] : 



— S hf _- h ee f 2z{l-z)Q—^- 



K (er) 
2vr 



(10) 
(11) 



where e 2 = z(l — z)Q 2 



m t 



The quark mass mf plays the role of a regulator when the photoproduction regime 



is reached. Namely, it prevents non-zero argument for the modified Bessel functions Ko t i(er) towards Q 2 —> 0. The 
electric charge of the quark of flavor / is given by eej. 

For vector mesons, the light-cone wavefunctions are not known in a systematic way and have to be determined in a 
phenomenological way. Here, we follow the approach proposed in Refs. [45|,|4(| and discussed in detail in [191 120I |54| . 
It is assumed that the spin and polarization of the vector mesons are the same as in the photon, which is considered as 
being predominantly a quark-antiquark state. Consequently, the wavefunctions for a transversely and longitudinally 
polarized vector meson are given by 



x +1 (T, z) = ±J ^ - 1 - f i 
Kh,\=±\\ > > V 2tt z(1-z) 1 



ie ±ier [z^ l ,±^ =F - (1 - z)5h, T 8- h ± ]d r + mf5 h ,±8j^ ± } 4> T {r,z). (12) 



and 



hh,\: 



=o( r > Z ) 



■ S h,-h 



M v + S 



V 2 



M v z(l - z) 



4>L{r,z), 



(13) 



where V 2 , = (l/r)d r + d 2 and My is the meson mass. As emphasized in (l9| . the longitudinally polarized wave function 
is slightly more complicated due to the fact that the coupling of the quarks to the meson is non-local, contrary to the 
photon case. 

Following Ref. (2lj we define the overlap function between the photon and vector meson wave functions by 



t>j' v (z,r;Q 2 ,M*) =^ [^(z,r;M 



(14) 



fhh 
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Using the above expressions for the photon and vector meson wave functions we get (See appendix of Ref. [HI): 

™2 _ V72 



<f>l" V (z,r, Q' 2 ) = e fx l^-N c 2QK Q (er) 



M v z{l- z)<j> L {r,z) + 



M v 



-<t>L(r,z) 



7*y 



(z, r, Q 2 ) = e f J ^N c ^f {m 2 f K (er)<t> T (r, z) - [z 2 + (1 - z^sK^drMr, *)} 



a e N t 



4tt 2tt 



(15) 
(16) 



where 



^l,t = N L:T exp 



m)R 2 
lz(l - z) 



m 2 R 2 



2z(l - z)r 2 
R 2 



(17) 



and the parameters R and Nl t & re constrained by unitarity of the wavefunction as well as by the electronic decay 
widths. The parameters used in our calculations are presented in Table[U Our choice is usually denoted in literature as 
boosted Gaussian (BG) wavefunctions and are a simplified version (5 = 1) of the wavefunctions proposed by Nemchik, 
Nikolaev, Predazzi and Zakharov in Refs. 45, 46]. We quote Refs. [l9l l20l. l2lj for more details and comparison with 
data for both photo and electroproduction of vector mesons in ep collisions at HERA. 
Using the above expressions for the vector meson wave functions, we obtain 



\XmA 



L,T 



d 2 bT A {b)\ \ I d 2 r I dz<S>lJ(z,r,Q 2 )o- dp exp{--<j dp T A (b)}\' 



and 



v coh (ll.T A -» VA) = Jd 2 bl^\Jd 2 rJ dz$l"%{z,r,Q 



2 [l-exp[--a dp T A (b)} 



(18) 



(19) 



In what follows we will estimate the total cross section for the exclusive vector meson production using as input in 
the Eqs. (TT5)) and (fTO)) the GBW and bCGC parameterizations for the dipole-proton cross section. 



III. RESULTS 



The project of an electron-ion collider is being elaborated by groups at BNL (eRHIC) 0] and at CERN (LHeC) 
[55l | . The main expectation is that the future eA collider should have a center of mass energy larger than 60 GeV 
and a high luminosity, such that measurements of inclusive observables at low- a; and diffractive processes could be 
improved substantially. Here we will analyze the energy and Q 2 dependence of the coherent and incoherent total cross 
section for J/ty and p vector meson production. We will consider two different nuclei (Ca and Pb ) in order to study 
the ^-dependence of our predictions. 

Moreover, we will also use as input in our calculations the linear limit of the GBW and bCGC dipole-proton cross 
section in order to determine the sensitivity of our results to the saturation effects in the proton. In the GBW model 
we take the linear limit by using instead of <r dp , the expansion of ([3]) when the dipolc is very small (r — > 0) and the 

exponent is small. In this case a dp = ctq M — . In the bCGC model the linear behavior is estimated by taking only 
the first line of Eq. © . 

Let us start our analysis considering nuclear photoproduction of vector mesons. In Figs. [1] and Owe present our 
predictions for the energy dependence of the coherent and incoherent total cross sections, respectively. As expected, 
the cross section grows with the energy and with the atomic number. We predict huge coherent cross sections, in 
reasonable agreement with those presented in [52], where a different approach was used to estimate the cross sections. 
At fixed atomic number, the coherent cross section is approximately a factor 70 (25) larger than the incoherent one for 
p {J/^f) photoproduction. Furthermore, at fixed energy, the coherent cross section increases with the atomic number 
faster than the incoherent one, in agreement with the conclusions from 26]. 

The photoproduction of p in nuclear collisions is a typical soft process, characterized by large dipole pair separations. 
At high energies, the saturation scale can assume large values and become the hard scale of the problem. The transition 
between these two regimes is determined by the energy dependence of the saturation scale, which is model dependent. 
From Figs. [1] and [2] we can see that for the p meson the full predictions of the GBW and bCGC models are very 
similar, in contrast with the linear ones. While the cross sections are strongly dependent on assumptions concerning 
saturation in the proton in the GBW case, they are very similar when we use the linear and full expressions of the 
bCGC model. This behavior is directly associated to the large difference between the linear and full expressions in the 




500 
400 
300 
200 
100 

800 °0 
5000 



1 1 1 
" J/ V P 


1 1 1 1 




bCGC 


Q 2 = 


— 




Ca 






Pb 






■ — ■ — Ca linear 






Pb linear 




' - ■ J . ■ - 


■■+-■■" 1"— 





200 400 600 800 



200 400 600 

W (GeV) 




800 



200 400 600 

W (GeV) 



800 



FIG. 1: Coherent photo-production cross sections for J/ip and p vector meson production with the GBW (left panels) and 
bCGC (right panels) models. 



GBW model at large dipole pair separation, which is smaller in the bCGC model. Concerning J /if? photoproduction, 
which is characterized by a hard scale, the mass of J/Sf?, the total cross sections increase faster with the energy than 
in the case of the p. Moreover, the GBW and bCGC predictions are very different, which is directly associated with 
the energy dependence of the saturation scale in these two models. Similarly to the p case, the linear and full GBW 
predictions are different, while are similar when we use the bCGC model. 

In Figs. [3]and[5]we present our predictions for the photon virtuality Q 2 dependence of the coherent and incoherent 
total cross sections, respectively. As expected, the total cross sections decrease with Q 2 . The Q 2 behavior predicted 
by the GBW and bCGC models are similar. The main difference is associated to the strong dependence on the 
saturation effects in the proton observed when we use the GBW model. Moreover, the Q 2 dependence of the J /if? 
cross section is smaller than the p one, which is directly associated to the hard scale present in the J /if? case. 

Finally, in Figs. [5] and [6] we present our predictions for the energy and Q 2 dependences of the ratio between the 
incoherent and coherent cross sections, respectively. The incoherent contribution is a small fraction of the coherent 
one and the ratio decreases with the energy. Moreover, the ratio is larger for small values of the atomic number. 
These conclusions are weakly model dependent. Concerning the Q 2 dependence, the ratio decreases at larger Q 2 , 
with a stronger dependence in the case of the p meson, as expected from the discussion above. It is important to 
emphasize that at large Q 2 the ratio is larger for J /if? than for p production. We can conclude that although measuring 
experimentally the intact recoil nucleus (coherent process) in a future electron ion collider may be very difficult, our 
results strongly suggest that this is the dominant process in vector meson production, specially when we focus on 
high energy eA reactions. This conclusion agrees with that obtained in 26]. 

So far we have been working only with two phenomenological models for the dipole-proton scattering amplitude, 
GBW and bCGC, and one might wonder what would be the results obtained with other existing models. In a 
previous paper we have performed a more systematic comparison between different dipole models (56j. Based on 
this comparison, we believe that using other scattering amplitudes would slightly change the normalization of our 




FIG. 2: Incoherent photoproduction cross sections for J/ip and p vector meson production with the GBW (left panels) and 
bCGC (right panels) models. 



predictions for the cross sections. However all the qualitative features would remain the same. In particular the low 
sensitivity to non-linear effects in the proton, observed in all the curves obtained with the bCGC model, would be 
the same for the other models. This can be traced back to the dependence of the saturation scale on the energy and 
dipole pair separation, which determine the speed with which the amplitude goes to one, which is very similar in the 
modern phenomenological models based on the CGC formalism. 



IV. CONCLUSIONS 



Although a large body of experimental data has been accumulated by the HERA collaborations HI and ZEUS 
over more than one decade in ep collisions, very little is known about vector meson production off nuclei at high 
energies. In this regime we expect a significant contribution from saturation physics. In this paper we have estimated 
the coherent and incoherent cross sections for the exclusive vector meson production considering the color dipole 
approach and phenomenological saturation models which describe the scarce data as well as the HERA data. Our 
results demonstrate that the coherent production of vector mesons is dominant, with a small contribution coming from 
incoherent processes. Probably in the future eA colliders the separation between coherent and incoherent processes 
will be difficult. However, in view of our results it might be worth trying. 




FIG. 3: Coherent electroproduction cross sections for J/ip and p with the GBW (left panels) and bCGC (right panels) models 
as a function of Q 2 at fixed center of mass energy (W = 100 GeV). 
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